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ABSTRACT 

We report on extragalactic sources detected in a 455 square-degree map of the southern sky made 
with data at a frequency of 148 GHz from the Atacama Cosmology Telescope 2008 observing season. 
We provide a catalog of 157 sources with flux densities spanning two orders of magnitude: from 
15mJy to 1500 mJy. Comparison to other catalogs shows that 98% of the ACT detections correspond 
to sources detected at lower radio frequencies. Three of the sources appear to be associated with 
the brightest cluster galaxies of low redshift X-ray selected galaxy clusters. Estimates of the radio 
to mm- wave spectral indices and differential counts of the sources further bolster the hypothesis that 
they are nearly all radio sources, and that their emission is not dominated by re-emission from warm 
dust. In a bright (> 50 mJy) 148 GHz-selected sample with complete cross-identifications from the 
Australia Telescope 20 GHz survey, we observe an average steepening of the spectra between 5, 20, 
and 148 GHz with median spectral indices of 05-20 = —0.07 ± 0.06, 020-148 = —0.39 ± 0.04, and 
a 5-i48 = —0.20 ± 0.03. When the measured spectral indices are taken into account, the 148 GHz 
differential source counts are consistent with previous measurements at 30 GHz in the context of a 
source count model dominated by radio sources. Extrapolating with an appropriately rcscalcd model 
for the radio source counts, the Poisson contribution to the spatial power spectrum from synchrotron- 
dominated sources with flux density less than 20 mJy is C Sync = (2.8 ± 0.3) x 10 _6 |xK 2 . 

Subject headings: surveys — radio continuum: galaxies — galaxies: active — cosmic background 
radiation 
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1. INTRODUCTION 

Large (> 100 square-degrees) millimeter- wave surveys 
are beginning to probe arcminute angular scales, corre- 
sponding to spatial frequencies t > 3000. At these small 
angular scales, the fluctuations in the extragalactic sky 
temperature become dominated by emission from galax- 
ies and the the rmal Sunyaev-Zeldovich (SZ ) effect from 
galaxy clusters (jSunvaev fc Zel'dovichllQTOt) , rather than 
primordial fluctuations in the Cosmic Microwave Back- 
ground (CMB). 

The predominant extragalactic point sources of emis- 
sion at 148 GHz (2.0 mm) are active galactic nuclei 
(AGN) and dusty, star-bursting galaxies. AGN de- 
tected at 148 GHz arc characterized by a synchrotron- 
dominated spectrum extending to lower radio frequen- 
cies. On the other hand, dusty, star-bursting galaxies at 
148 GHz display a grey body spectrum increasing with 
frequency into the submillimcter. The source of 148 GHz 
flux from AGN is synchrotron emission concentrated near 
the central accreting super massive black hole while for 
dusty star-bursting galaxies the millimeter flux is sourced 
by re-emission from dust that is heated primarily by 
prodigious star formation. While the millimeter emission 
from the majority of dusty star-bursting galaxies is be- 
low the nominal flux density limits of curre nt large-scale 
survey s, a recent 150 and 220 GHz study bv lVieira et al.l 
(|2010h using the South Pole Telescope identifies a sub- 
population of these sources with anomalously high fluxes 
which likely belong to a rare, lensed population of high- 
redsh ift dusty galaxies (jNegrello et all 120071 : iLima et all 
2010). It follows that current and future wide-area mil- 
limeter surveys will identify important sub-populations 
of core-exposed, radio-loud AGN and lensed dusty galax- 
ies. With the members of these sub-populations identi- 
fied, detailed follow-up studies will help us better un- 
derstand these high-energy states of galaxy formation as 
well as their important role in providing energy feedback 
to their environments. 

Because source emission is a significant contributor to 
the overall sky brightness at small scales, the charac- 
terization of extragalactic sources is essential for inter- 
preting the primary CMB anisotropics and the SZ signal 
from galaxy clusters. Measurements of the primary CMB 
power spectrum at high spatial frequencies (I > 2000) 
will constrain the form of the inflationary potential (e.g., 
the spectral index n s of primordial fluctuations). Such 
measurements will require information about the spec- 
tral and spatial distribution of millimeter sources in or- 
der to separate foregrounds from the primordial signal. 
At still smaller scales {I > 3000), studies of the CMB 
spectrum attempt to constrain fluctuations in the mat- 
ter density field from the contribut i on of the SZ to th e 
power spectrum dLueker et alll2010 l : iFowler et al.ll2010f) . 
For these studies an in-depth understanding of the point 
source populations is even more critical for separating 
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the power spectrum of the SZ from that of sources. Fur- 
thermore, an understanding of the energy feedback from 
AGN and star-formation to the cluster environment will 
be im portant for constrain i ng the form of the SZ spec- 
trum (JBattaglia et alJl2010HShaw et alJl2010HTrac et all 
201 Q). Finally, SZ surveys attempting to measure Hm, 
as, and dark energy through the evolution of the clus- 
ter mass function will likewise need to consider the spec- 
tral behavior and cluster occupation numbers for sources 
in order to avoid systemat ically biasing mass estimates 
based on SZ flux density dLin et all 12009 : Sehg al et all 
l2lMlVanderlinde et al.H2010D . 

The Atacama Cosmology Telescope (ACT) is a 
millimeter-wave observatory which will ultimately sur- 
vey thousands of square degrees of sky at arcminute res- 
olution with milli-Jansky sensitivity to sources. ACT is 
located at 5200 m in the Atacama Desert in the Andes 
of northern ChilcQ This high desert site in the trop- 
ics was chosen for its excellent atmospheric transparency 
and its access to northern and southern celestial lati- 
tudes. ACT observes simultaneously in bands centered 
at 148GHz (2.0mm), 218 GHz (1.4 mm), and 277 GHz 
(1.1 mm), each band having a dedicated 1024-clcmcnt 
array of bolometric transition edge sensors. As of mid- 
2010, ACT has completed three seasons of observations: 
2007, 2008, and 2009, and the 2010 season is underway. 
In each season, ACT has conducted two surveys: a 9°- 
wide survey centered at declination —53.5° and a 5°-wide 
stripe centered on the celestial equator. 

In this paper we report on extragalactic sources in the 
ACT 2008 148 GHz dataset. This is the first report 
devoted to ACT source science a nd complements th e 
148 GHz power spectrum study in I Fowler et all (|2010[) . 
In what follows we give an overview of the observations 
and data reduction (Section [2]) , describe the source cat- 
alog (Section |3] and the Appendix), and discuss implica- 
tions of the study including constraints on source models 
(Section |U). 

2. OBSERVATIONS AND DATA REDUCTION 

The data used for this analysis were collected by ACT 
at 148 GHz during its second observing season in 2008. 
This section gives an overview of the survey observations 
and the reduction of the raw data to a map as well as 
a detailed description of the source extraction. For a 
more thorough introduction to the ACT facility, obser- 
va tions, and da t a red u ction pipeline, we r efer the reader 
to lFowler et~all (|2010D . lSwetz et all (120101 ) and references 
therein^ 

2.1. Observations 

The 2008 southern observations were carried out over 
a survey area 9° wide, centered on declination —53.5°, 
and extending from right ascension 19 h to 24 h and 
00 h to 07 h 30 m . The subset of these data used in the 
present analysis lies between right ascensions 00 h 12 m 
and 07 h 08 m and declinations -56°11' and -49°00' (455 
square-degrees). The area was chosen t o encompass 
the data used for power spectrum work in IFowler et all 

1 The ACT Site is at 22.9586° south latitude, 67.7875° west 
longitude. 

2 ACT C ollaboration papers arc archived at 
http://www.physics.princeton.edu/act/ 
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Fig. 1. — Sensitivity map with detections. The subset of the ACT 2008 148 GHz dataset considered for this study lies between right 
ascension 00 h 12 m and 07 h 08 m and declination — 56°11' and — 49°00' (455 square-degrees). The gray-scale encodes the rms of the map in 
mjy. The deepest data correspond to an exposure time of 23.5 minutes per square-arcminute and a 1<t sensitivity of 2.5 mjy. White circles 
mark the locations of ACT sources. The diameter of each circle is proportional to the log of the associated source flux density. Towards the 
edge of the map, the noise properties display local variation. For this reason, detections with flux density values below 50 mjy have been 
discarded in regions where the rms exceeds 4.6 mjy, corresponding to less than 7 minutes of integration. This, toget her w ith an exclusion 
of all detections below 5.25 cr, accounts for the relative dearth of detections in areas of shallow coverage (See Scction l2.ll ). 



(|2010( ) and represents a large fraction of the most deeply 
covered regions from the 2008 148 GHz dataset. Figure 
Q] shows the area of the sky used and associated point 
source sensitivies. Typical white noise levels in the map 
are 30-50 |j.K-arcminute, tending to higher values to- 
wards the map boundaries. As described in Section 12.41 
this white noise level, when match filtered with the ACT 
beam, results in typical sensitivities to point source flux 
densities from 2.5 to 5 mjy. 

The 2008 ACT observing season extended from mid- 
August to the final week of December. Observations 
took place during nighttime hours: from roughly 20:00 
to 06:00 local time. Of the total observing time, approx- 
imately 85% was devoted to the southern region. ACT 
observed by scanning at a constant elevation of 50.5° 
while the survey region drifted through the scan with the 
rotation of the Earth. During the first half of a night, 
ACT scanned at azimuth 150°, targeting a rising section 
of the survey area. During the second half of the night, 
ACT scanned the same section setting on the other side 
of the south celestial pole at azimuth 210°. The rising 
scans cross the survey region from southwest to northeast 
and back (in equatorial coordinates), while the setting 
scans cross the survey region from southeast to north- 
west. Together, the rising and setting scans cross- link 
each point on the sky with all adjacent points. The re- 
sulting cross-linked temperature data in principle contain 
all information necessary to recover an unbiased, low- 
noise map of the millimeter sky. In addition to survey 
observations, ACT also executed regular observations of 
Uranus and Saturn during 2008 to provide gain calibra- 
tion, beam profiles, and pointing. 

With the telescope scan strategy described above, any 
given location in the survey area would be observed over 
a period of approximately two months during a season. 
Therefore, source flux densities reported here are the av- 
erage flux density over a two month period. This is an 
important point as the vast majority of sources presented 
in th is paper are AGNs which are known to be variable 
(e.g.. lKesteven et al.iri977HValtaoia et al.lll992[ ). 



2.2. Reduction to Maps 

The raw 148 GHz data consist of 1024 time-ordered 
data streams, one per element of the detector array. Ap- 
proximately 25% of the data are rejected on the basis 
of telescope operation and weather. After further cuts 
based on individual detector performance, the data from 
680 148 GHz detectors over 850 hours (~ 3200 GB) are 
retained from the 2008 southern survey. 

Pointing reconstruction is accomplished in two steps. 
First, the relative detector pointings arc established with 
1.2" certainty through observations of Saturn. Second, 
absolute detector array pointings for our two southern 
survey configurations (rising and setting: 150° and 210° 
azimuth, 50.5° elevation) are established with 3.5" pre- 
cision through an iterative process in which the absolute 
pointing is adjusted based on offsets of ACT-observed 
radio source locations with respect to source locations 
take n from the Australia Telescope 20 GHz (AT20G) sur- 
vey (jMurphv et al.ll2010f ). 

Nightly calibrations of the detectors' responsivity 
(power-to-current conversion) are based on load curves 
taken at the start of each night. Stability of this cali- 
bration through the night is monitored using small steps 
in the detector bias voltages and established at the few 
percent level. Relative detector flux density calibrations 
are based on normalizing the detector responses to the 
beam-filling atmospheric signal. The resulting relative 
calibration is shown to be constant through the season 
at the few percent level. The final brightness tempera- 
ture calibration is based on ACT observations of Uranus 
througho ut the season and t he WMAP7 Uranus tem- 
perature (|Weiland et al.ll2010l ) extrapolated to the ACT 
14 8 GHz band. The calibration is more fully described 
in lFowler et al.l (I2010D . The overall calibration is certain 
to 6% rms, a number dominated by systematic uncer- 
tainties in extrapolating the temperature of Uranus to 
148 GHz from WMAP frequencies. 

The final step in the data reduction is map-making. 
An iterative preconditioned conjugate gradient solver is 
used to recover the maximum likelihood (ML) maps. The 
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Fig. 2. — Model ACT auto-power spe ctra used for the matched filter. Left. One-dimensional spectra decomposed by component. The 
CMB spectrum is taken from WMAP5 fNolta ct al. 2009). The undetected s ources and the Sunyaev-Zel'dovich effect from galaxy clusters 
are ID templates fit to the high-£ f 48 GHz spectrum in [Fowler °t al. (20l3). The fD noise spectrum is obtained by radially-binning the 
average of 2D spectra from ACT f48GHz jackknifc maps. The spectra are dominated at low-£ by the CMB and atmospheric noise and at 
high-£ by white detector and photon shot noise. Convolution of all celestial components by the ACT beam results in attenuation of the 
corresponding spectra at high-^. Right. The model two-dimensional auto power spectrum. The spectrum includes noise, CMB, undetected 
sources and the Sunyaev-Zel'dovich effect from galaxy clusters. The latter three components were obtained from the ID models. The 
ACT noise is isotropic except for extra noise in scan parallel (from residual 1/f) and scan perpendicular (from detector row correlations) 
directions. For each rising and setting scan direction, two orthogonal bands of excess power are centered on the origin rotated 60 with 
respect to one another. To minimize contamination, it is important to deweight data in these diagonal modes through the two-dimensional 
matched filter. 



algorithm solves simultaneously for the millimeter sky as 
well as correlated noise (e.g., a common mode from atmo- 
spheric emission). The map projection used is cylindrical 
equal area with a standard latitude of —53.5° and 0.5' 
square pixels. For more details on the mappi ng and other 
reduction steps, refer to lFowler et al.l (j2010[ ). 

2.3. Data Modeling 

The next step toward a source catalog is the construc- 
tion of a filter which optimizes the signal-to-noise ratio 
(S/N) of sources in the 148 GHz mapfj In order to con- 
struct such a filter it is necessary to obtain best estimates 
of the power spectra of the different components, signal 
and noise, which contribute to the ACT data. We model 
the temperature T at position x as the sum over sources 
plus other components in the map: 



T{x) = ^Tib{i 



T, 



other 



(x) 



(1) 



where Tj is the peak amplitude of the i source, b is the 
ACT 148 GHz beam function n ormalized to unit a mpli- 
tude and taken to be isotropic (jHincks et al.ll2009t ). and 
T'other includes contributions from the primary CMB, un- 
detected point sources, SZ from clusters and noise from 
the detectors and atmosphere. 

Figure [2] shows graphs of the power spectra of the 
contributors to T ot her- For the primary C MB compo- 
nent w e use the WMAP5 best fit model from lNolta et al.l 
(|2009t ). The SZ and undetected source components are 
taken f rom models fit to th e ACT 148 GHz power spec- 
trum inlFowler et al.l (120101) . The spatial power spectrum 
from iFowler et al.l ( 20101 ) was computed after masking 



See Section 12.41 for the definition of S/N used in this work. 



sources from a preliminary version of the catalog pre- 
sented in this paper. As such the derived contribution 
to T thcr from undetected sources should be a suitable 
estimate for the matched filter in this work. Specifically, 
the source model is a Poisson spectrum normalized to 
e(£ + l)C p /2ir = 11.2 u.K 2 at I = 3000, and the SZ model 
is a combined the r mal a nd kinetic SZ model (os = 0.8) 
from lSehgal et all (J2010T) normalized by a best- fit factor 
of 0.63. These celestial models (CMB, sources, and SZ) 
are all convolved with the ACT 148 GHz beam function. 
The celestial models are natively one-dimensional func- 
tions of multipole £ (as shown in the left graph in Fig- 
ure [2]) with corresponding azimuthally symmetric two- 
dimensional representations (as shown in the right graph 
in Figure [2]). The noise term, on the other hand, is na- 
tively a two-dimensional power spectrum estimated by 
the average of power spectra from difference maps. The 
difference maps were constructed by subtracting a map 
made from one half of the data (with a random selection 
of observing days contributing) from a map made from 
the other half (using the remaining observing days). The 
two-dimensional noise power spectrum is binned radially 
to obtain the one-dimensional representation shown by 
the line of connected dots in the left graph of Figure [2] 

From Figure [2] it is clear that the primary contributions 
to Tother come from the CMB and atmosphere at large 
scales (£ < 2500) while the white detector and photon 
shot noise dominates the power at small scales. In the 
spectral trough around i — 2500 — 3500 the detectable 
sources will have their greatest signal-to-noise ratio. 

An important feature of the two-dimensional power 
spectrum in the right graph of Figure[2]is the anisotropic 
nature of the noise term which corresponds to striping in 
the map. The stripes are a result of large scale drifts in 
atmospheric emission along the scan directions as well 
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Fig. 3. — 148 GHz map. The submap shown above is a sample 64 square-degrees of the survey field. The data have been matched 
filtered such that the gray-scale is in units of flux density (m Jy ) . The inset shows the flux density distribution across the data weighted by 
the max- normalized square-root of the number of data per pixel yj N \, s (x) / ^/JV t,s !max . Thus the distribution represents the data in the 
deepest part of the map (although it uses weighted data from all regions) . The data distribution is shown as a grey histogram on which 
is plotted a dashed Gaussian distribution with standard deviation 2.5 mjy. The positive non-Gaussian tail may, in p art, be attributed 
to sou rces and the negative tail to ringing from the filter about sources as well as SZ (e.g., ACT-CL J0438-5419 from Mcnantca u et al.l 
1120101 )). Several sources are apparent as white points surrounded by darker rings from the filter. The white contour marks the transition 
at the edge of the map where the rms exceeds 4.6 mjy, corresponding to less than 7 minutes of integration per arcminute. In this region 
we have excluded detections below 50 mjy due to contamination from local noise. The source above the contour at <5 ss —50° has a flux 
density of 150 mjy and is included in the catalog. 



as from correlations among rows of ACT detectors per- 
pendicular to the scan directions. Accounting for this 
anisotropy when filtering the map is important for ex- 
tracting an uncontaminated sample of sources. To prop- 
erly down-weight these noisy diagonal modes, we adopt a 
two dimensional noise covariance in the matched filtering 
technique described in Section l2~4l 

The final significant noise term that is not captured 
in difference maps is consistent with scan synchronous 
noise. The noise is likely due to instabilities induced by 
acceleration at scan turn-arounds. The noise manifests 
itself as horizontal striping in the maps and as excess 
power isolated to a vertical strip —100 < £ x < 100 in 
Fourier space. For simplicity, we set the noise power to 
infinity in this section of Fourier space such that modes 
contaminated by this noise will be completely down- 
weighted by the filter as described in the next section. 



2.4. Detection and Catalog Construction 

To optimize the S/N of the sources with respect 
to the background, we use a matched filter. This 
approach has been pro posed and used in previous 
work t o find both sources (iTegmark fc de Oliveira-Costal 
[T998L iWright et all I2FM IVieira et al.l I2010D and the 
Sunyaev-Zel'dovich effect fro m clusters ( Melin et al.l 
120061 : iStaniszewski et aLll2009D . For completeness of pre- 
sentation, here we rederive the form of the filter and its 
basic properties. Without loss of generality, we consider 
a field centered on a source such that 



T{x)=T Q b(x)+T other {x). 



(2) 



We apply a filter $(fc) in Fourier space such that the 
temperature at the center of filtered field is 



Tmt(0) - T / $(fc)6(fc)dfc 



J $(fe)T c 



other (fc)dfc (3) 



6 



where b and T ot hcr are the Fourier transform of the beam 
and noise temperature field, respectively. The S/N of the 
central source is 



S/N 



T / <S>(k)b{k)dk 



J$(fe)T other (fe)dfe 

A filter that maximizes this S/N is 

b*(k) | f other (fe) |- ; 



$(fc) = 



/&*(fc')|T othor (fc')h 2 &(fc')dfc' 



(4) 



(5) 



where the normalization has been chosen to produce an 
unbiased estimate of the amplitude of the source To at 
x = 0. The noise variance of the filtered data is 



*(fc)T ot hcr(fe) \ Z dfc 

b*(k) | T othcr (fc) |- 2 b(k)dk 



(6) 



In practice, the noise variance of the filtered map is 
obtained from the filtered map itself after masking the 
brightest six sources (S/N > 50). These sources increase 
the rms of the maps by 10%. The remaining sources 
contribute approximately 1% to the rms. This slightly 
more conservative estimate of the variance agrees with 
the estimate from Equation [5] in which T'other has been 
constructed as described in Section [2751 In what follows, 
the S/N at a location in the filtered map is defined as the 
temperature at that location divided by the square-root 
of this variance. 

Before applying the global matched filter from Equa- 
tion [5j we multiply the map, pixel-wise, by the square- 
root of the number of observations per pixel normalized 
by the observations per pixel in the deepest part of the 
map, -\/A' r o bs(a0/-ZVobs,max- This is equivalent to weight- 
ing the data by the inverse of the estimated white noise 
rms shown in Figure [Hand accounts for local variation in 
the white noise amplitude. Furthermore, the map is ta- 
pered to zero in a 10' boundary region around the edge 
of the map to mitigate the artifacts arising from data 
aperiodicity when filtering. This windowed data is ex- 
cluded from the final analysis, reducing the usable sky 
area from the total 455 square-degrees to 443 square- 
degrees. Next, because the ringing of the filter around 
the brightest sources can cause false detections, we iden- 
tify and mask the six most significant ( > 50ct ) sources 
before applying the matched filter to the map. These 
brightest sources are treated and included in the catalog 
in the same manner as the fainter sources with the ex- 
ception they are recovered through an initial run of the 
pipeline with the S/N lower limit increased. An impor- 
tant final caveat: in constructing the noise term T ot hcr, 
the component models need to be tapered and weighted 
in the same fashion as the data for the matched filter 
formalism to hold. This is particularly true for compo- 
nents with red spatial spectra, such as the atmosphere 
and CMB, because aliasing due to a particular window- 
ing scheme can significantly alter the spectrum. Figure 
[3] shows a sample 64 square-degrees of the filtered map. 
For this reason the catalog includes only detections with 
S/N 5.25 and above. See Section [3~3l for a discussion of 
purity determination. 



Localized, non-white noise in the map requires that 
we take further measures beyond the global matched fil- 
ter solution outlined above. First, local large-scale at- 
mospheric noise requires us to add a low-^ taper to the 
term | T thcr(fe) | -2 (Equation [5]) which rises from zero 
at t = to one at £ = 1200 as sin 5 (7r^/2400). This 
filter removes the local atmospheric noise while down- 
weighting only a small fraction (~1%) of the data con- 
taining source power in Fourier space. Second, in areas of 
the map which are particularly shallow, uneven coverage 
leads to excess striping. In these parts of the map the 
noise model described in Section 12.31 is invalid and non- 
white noise remains even after the filter is applied. For 
this reason, we exclude sources from the catalog which 
are detected with flux densities below 50 niJy in areas 
of the map with exposure times less than 7 minutes per 
square-arcminute. This exposure time is approximately 
one-third the exposure time in the deepest areas of the 
map and corresponds to 4.6 mJy rms. The cut level of 
50 mJy at 7 minutes per square-arcminute was chosen to 
broadly eliminate contamination observed in simulations. 
Future studies will make use of local noise estimation to 
avoid such an exclusion. The sample submap in Figure[3] 
demarcates the region at the edge of the map in which we 
exclude 50 mJy detections with a white contour. With 
this exclusion, the area used for sources with flux density 
below 50 mJy is 366 square-degrees. 

The final step in the catalog generation is to derive the 
flux densities associated with the detections. Given the 
form of the filter in Equation^ the source-centered value 
of the filtered map, multiplied by the solid angle of the 
beam profile, is the source flux density. It is this value, 
rescaled by the inverse of the square-root of the num- 
ber of observations per pixel normalized by the number 
of observations per pixel in the deepest part of the map 
(\/ N max /N bs,max(x)), that we record as the raw flux 
density estimate for a detection. It follows that an error 
in source location results in an error in flux estimation. 
Such an error arises due to finite pixel size: a detection 
rarely falls in the center of a pixel. Source location er- 
ror due to finite pixel size causes a systematic negative 
bias and increased scatter. The 0.5' ACT map pixels 
cause a 10% negative bias with comparable scatter in 
the flux density estimate. To remove this bias and scat- 
ter, we zero - pad t he filtered data in Fourier space (e.g., 
I Press et all (|1992f) ) such that the pixel spacing in map 
space is decreased by a factor of sixteen. This Fourier 
interpolation of the filtered data onto 0.03125' pixels is a 
convenient way to better locate the peak of source emis- 
sion and therefore mitigate systematic errors in position 
and flux density estimation to less than 1%. An impor- 
tant caveat to this technique is that the pixel window 
function must be taken into account when reducing the 
pixel size. As discussed further in Section |3.2.3[ simula- 
tions show that flux densities thus derived are unbiased 
at the sub-percent level. 

3. CATALOG 

The ca talog of 157 ACT extragalactic sources is given 
in Table IA1I of the Appendix. The catalog provides 
the IAU name, celestial coordinates (J2000) as well as 
S/N and 148 GHz flux density estimation of each ACT- 
detected source. Raw flux densities arc estimated di- 
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Fig. 4. — Astrometric Accuracy. The small filled circles are the 
positional offsets of ACT sources with S/N > 20 from counter- 
parts in AT20G. The error bars show the rms in right ascension 
(3.5") and declination (3.3") and are centered on the mean of the 
distribution: —0.3" ± 0.7" in right ascension and 0.0" ± 0.6" in 
declination. 



rectly from the map as described in Section 12.41 De- 
boosted flux densities (Section r3.2.2p are also given. The 
seventh column gives the exposure time in minutes per 
square-arcminute at the location of the source. Finally, 
we provide the ID of AT20G sources collocated within 
30" of an ACT source. The following subsections pro- 
vide the necessary information for interpreting the cata- 
log data. 

3.1. Astrometric Accuracy 

The positions of the twenty-seven ACT sources with 
S/N > 20 were compared to positions of associated 
sources in the AT20G catalog. The AT20G pointing 
is checked against VLBI measurements of International 
Celestial Reference Frame calibrators an d is shown to 
be ac curate to better than one arcsecond ([Murphy et al.l 
120101) . Figure 2] shows the offsets of the twenty-seven 
ACT source locations with respect to associated AT20G 
locations. The mean of the offsets is —0.3" ± 0.7" in 
right ascension and 0.0" ± 0.6" in declination. The rms 
of offsets is 3.5" in right ascension and 3.3" in declina- 
tion. One ACT source with S/N > 20 is cross-identified 
with an AT20G source 18" away. This source is ACT-S 
J042906-534943 and is flagged as extended in the AT20G 
catalog. Furthermore, there are multiple SUMSS sources 
located within the search radius suggesting a nontrivial 
geometry. Therefore, this source represents an outlier 
in this analysis and is not used. For twenty-four cross- 
identified sources with S/N < 8, the ACT location rms 
with respect to AT20G position was inflated by the effect 
of the noise: 8.8" in right ascension and 6.2" in declina- 
tion. 

We have identified a systematic shift in pointing 
throughout the season likely associated with movement 
of the telescope. A significant fraction of the scatter in 
FigurcS]is attributable to this systematic effect. Because 
of this effect, the distribution of pointing errors in Figure 
[4] docs not appear Gaussian. We plan to eliminate this 
effect in a future study. 
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Fig. 5. — Source Profile Fitting. The radially binned and nor- 
malized brightness profile observed for ACT-S J021046-510100 is 
represented by the points with error bars. The solid line running 
throug h the points is the measured ACT beam from lHincks et all 
H2OQ90 convolved with a Gaussian with a = 3" plus an additive 
background term. The Gaussian convolution is intended to repre- 
sent spreading of the beam by random pointing error. The value 
of cr = 3" describes the best-fit to a Gaussian-convolved beam plus 
background term. Note that 3" is consistent with pointing uncer- 
tainty from Figure [4] The dashed line shows the model with back- 
ground and beam amplitude fit to the data, but no Gaussian con- 
volution. Without the Gaussian convolution, the background term 
increases to fit the data at small angles and forces the model high 
at large angles. The data clearly prefer the Gaussian-convolved 
model. Similar fits to the twenty brightest sources in the ACT 
data suggest bias in the reported flux density level due to misesti- 
mation of the source profile is at the sub-percent level. The error 
bars shown on the radially binned points are correlated. 



3.2. Flux Density Recovery 



As discussed in Section [2721 the overall calibration has 
an uncertainty of 6%, dominated by systematic errors 
in the temperature of Uranus. In addition to this un- 
certainty, errors in the flux estimation may arise due to 
errors in the assumed source profile, to flux boosting of 
lower significance candidates, and to a failure of the map- 
maker to converge. In this section we describe tests of 
these potential sources of flux density error as well as 
an end-to-end check of the match-filter recovery and de- 
boosting through simulations. 

3.2.1. Beam Profile 

Flux density recovery is a function of the form of the 
source profile b assumed in the filter (See Equations [5] and 
[3). For the filter's source profile b we a dopt the ACT 
148 GHz beam from IrTincks et al.l (|2009i ). Deviation of 
the actual source profile from the ACT 148 GHz beam 
will result in a biased estimate of the flux density as 
determined from the filtered map. 

To search for a difference between the actual and as- 
sumed profiles, we examine the profiles of the twenty 
most significant sources in the ACT map. We fit the 
data with a background term plus an altered ACT beam 
(respectively broadened or squeezed by the convolution 
or deconvolution with a Gaussian). The Gaussian convo- 
lution wc employed models the broadening of the beam 
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Fig. 6. — Flux density dcboosting. The Gaussian flux density 
distribution P(S m \ So) (dashed) derived directly from the maps 
is multiplied by the prior distribution P(So) (dotted) of the intrin- 
sic source flux to obtain the deboosted posterior flux distribution 
(solid). Since the mean of the maps is set to zero through filtering, 
P(S'o) can assume negative values. The peak in the distribution 
shifts to lower density flux (a <10% effect for 20mjy sources). 
For the S/N=5.6 source shown here, there is little volume in the 
posterior P(S m | So) at negative flux. 



resulting from the random pointing error. The average 
best-fitted profile is the beam from iHincks et all (J2009I ) 
convolved with a Gaussian with a = +3" ± 1". This 
smearing is consistent with the pointing rms in Figure 01 
The fit is performed out to a radius of 3'. Figure [5] shows 
the fit for the brightest source. The effect of this convo- 
lution on the 84" full width at half maximum (FWHM) 
of the ACT beam is at the sub-percent level: 0.05". Flux 
density misestimation will scale roughly as the ratio of 
the actual profile solid angle to the assumed profile solid 
angle. In the present case, the misestimation is below 
1%. 

As an additional check on the reported flux densities, 
the matched filter based flux density estimates are com- 
pared to l'-diameter aperture flux densities. The aper- 
ture flux density is defined as the background subtracted 
flux within a circle (the aperture) centered on the source. 
The aperture flux density is estimated from the unfiltcred 
ACT maps and as such is complementary to the flux den- 
sity estimates from the filtered maps. As the source pro- 
file extends well beyond a 1' diameter, the aperture flux 
density estimates are expected to be biased low with re- 
spect to true flux density. For the thirty brightest sources 
from the ACT data, the average ratio of aperture flux 
density to matched filter flux density is 0.610 ± 0.003. In 
simulations which featured sources with the ACT beam 
shape, the same ratio is 0.607 ± 0.004. If the source 
profiles in the data were different from the ACT beam 
shape, then this ratio from the simulation would differ 
from the data. The agreement between data and sim- 
ulation further bolsters the claim that flux density bias 
due to source profile misestimation is at the sub-percent 
level. 

3.2.2. Deboosting 

In a source population for which the counts are a 
steeply falling function of flux, a source's measured flux 
density S m in the presence of noise is likely to be an over- 



estimate of its intrinsic value Sq. The overestimation is 
most pronounced for sources at low S/N. The process of 
deboosting accounts for this overestimate by construct- 
ing the posterior flux distribution based on prior knowl- 
edge of the source population. Given the relatively high 
flux density and significance of the detections in this 
work ; we adopt the strai ghtforward Bayesian approach 
from lCoppin et all pOOl : 



P(S \S m )cxP(S m \S )P(S ) 



(7) 



where the probability P(S m \ Sq) of measuring S m given 
So is taken to be normal with mean Sq and variance 
derived from the S/N. 

The prior probability of flux Sq in a pixel, P(Sq), 
is computed by generating simulations of the (filtered) 
intrinsic source flux distribution per pixel. Individ- 
ual, filtered source profiles, To J $(fc)6(fc)dfc, arc added 
to a blank survey map at randomly chosen locations. 
The numbers and associated amplitudes, To, of the 
sources are cho sen in accordance wit h infrared and ra- 
dio counts from lToffolatti et al.l (|1998[ ) rescaled to fit the 
counts from this study (See Section 14.41 ). The proce- 
dure for computing the prior probability based on source 
counts fit to ou r data involved an initial rescaling of the 
iToffolatti et al.l (|1998l) data model to counts computed 
with raw flux densities. Given deboosted flux densi- 
ties based on this initial prior, counts were then recom- 
puted, the prior was re-estimated and and used to obtain 
new deboosted flux densities. This process was iterated 
until corrections to the counts and the deboosted flux 
densities became negligible. Only one iteration was re- 
quired for the relatively small level of deboosting used 
in this initial study. Furthermore, we cut off the radio 
counts at 150 mJy, reflecting the fact that the brightest 
six sources (S > 150 mJy) are detected and subtracted 
before constructing the rest of the catalog. A function 
Pi(So) is then generated by binning the fluxes associated 
with map (indexed here by i) pixels in 0.5 mJy bins. 
The final distribution P(Sq) was then computed as the 
average of Pi(So) from ten-thousand independent sim- 
ulations. The deboosting algorithm is illustrated for a 
single source in Figure [5] where the dashed Gaussian rep- 
resents P(S m | Sq) and the dotted profile peaking just 
below zero flux is the prior probability P(Sq). The pos- 
terior probability P(So\S m ) is the solid line. 

The deboosted flux Sdb reported in the ACT catalog 
for sources below 50 mJy is the median of the associated 
P(So\S m ), and the reported asymmetric errors enclose 
the 68% confidence interval. The abrupt 150 mJy cutoff 
imposed on the radio counts in combination with finite 
pixel size effects the smoothness of the prior estimate 
P(Sq) at higher fluxes. To mitigate this effect, the sim- 
ulated maps used to construct the prior feature a pixel 
size of half that of the data (0.25'). Furthermore, de- 
boosted fluxes are only provided for sources with flux 
below 50 mJy where the computed prior is smooth: for 
sources above 50 mJy, we simply report the center and 
68% confidence level of P(S m \S ). 

The prior probability P(Sq) in Figure [6] is broader and, 
for the range of flux densities plotted, more symmetric 
than analogous di stributions d erived in previous work 
(e.g., Figure 6 of iScott et all (12001 ). The difference 
arises because our simulations include the radio popu- 



lation whose source counts are much shallower than the 
infrared populations. The bright radio sources, in addi- 
tion to having a bright positive tail, produce significant 
negative ringing when filtered. A more familiar form of 
P{Sq) would be generated if we were to rec over fluxes us- 
ing the CLEAN technique (|Hogbomlll974[ ) as described 
in, e.g.. lVieiraet all (J2010D . 

In future multi-band work that considers detections 
of lower significance, we will employ important ex- 
tensions of the s imple deboosting used her e, e.g., 
lAustermann et al.1 (|2010Q : ICrawford et all (|2009| ). 



TABLE 1 
Number Counts, Purity, and Completeness. 3 




S (mjy) 

Fig. 7. — Simulation Flux Recovery. The match- filter derived 
and deboosted flux density estimate S^b plotted versus the intrinsic 
flux density. The dashed gray lines show the ±15 mjy limits about 
the solid black one-to-one line. For So < 20mjy, the sample is 
incomplete due to the S/N > 5.25 selection. This is manifest in 
the graph by the apparent skew of the population above the one- 
to-one line for low flux densities. 



3.2.3. Simulations and Flux Recovery 

Having established that the source profile used in the 
matched filter is a good approximation to that in the 
data (Scction l3.2.II) . it remains to test the flux estimation 
and deboosting through simulation. Source flux densities 
were derived from maps with celestial components and 
noise modeled as described in Section l2"T3"l with two excep- 
tions. First, the SZ component was excluded to prevent 
confusion due to cluster-galaxy correlations. This effect 
caused negative outliers in derived flux density due to 
brightness cancellation by the cluster decrement. Sec- 
ond, we excluded several (> 10) peculiar instances of 
superposed, very bright, (30 mjy) infrared sources. This 
apparent clustering of ultra-bright dusty sources caused 
positive outliers in derived flux density. 

With the exceptions noted above, the simulated 
sources were detected with the same blind algorithm 
that was described for source extraction from the data. 
Similarly, the deboosted flux densities of the simulated 
sources were derived. Figure [7] shows the result. At the 
percent level, the derived deboosted flux density Sdb is a 
consistent estimate of the intrinsic flux density So with 
a best fit of 
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-0.09 ± 0.29mJy + 1.008 ± 0.005 x S (8) 



where the fit (reduced x 2 — 1.06 for 451 degrees of free- 
dom) was performed only for sources with flux density 



a See Figure lTul for a graph of purity/complcteness-corrected 
differential source counts. 

b Fo r a discussion of flux density dependent area, see Section 
12^1 

greater than 20mJy to restrict the analysis to a com- 
plete (and thus symmetric) distribution. The errors (68% 
c.l.) are derived from one thousand bootstrap samplings. 
Furthermore, the model Sdb = So fits the 451 sources 
with flux densities greater than 20mJy with a reduced 
X 2 = 1.07. This one-to-one model fits the entire popu- 
lation (down to the S/N=5.25) with reduced x 2 = l-19 
over 651 degrees of freedom. This improbable statis- 
tic likely results from a combination of unaccounted- 
for as trophysics (e.g., source clustering) in iSehgal et al.1 
(|201CO . errors in our deboosting, and an underestimate 
of flux density error. 

3.2.4. Convergence 

As described in Section [2] maximum-likelihood maps 
generated from the cross-linked ACT 148 GHz data are 
unbiased for modes corresponding to multipoles in excess 
of a few hundred (i.e., exactly those used for source flux 
estimation). We checked this claim through signal-only 
simulations and found that the flux density estimates of 
the sources converged to within 1% of their simulation 
values before the twenty-fifth iteration of the precondi- 
tioned conjugate gradient map-solver (See iFowler et al.1 
<|2010f) for a fuller description of the solver.). This study 
uses a map from the hundredth iteration of the ACT 
maps. Therefore, we conclude that the flux densities are 
not biased by failure of the mapping algorithm conver- 
gence. 

3.3. Purity and Completeness 

The number of false detections in our catalog of 157 
sources is estimated by running the detection algorithm 
on an inverted (negative temperature) map in which 
the SZ decrements from all ACT-detccted and optically- 
confirmed clusters have been masked. A full description 
of the ACT SZ cl uster population and op tical follow- 
up can be found in iMenanteau et "ail (|2010D . With this 
approach, three spurious detections are found, giving a 
purity of 98% for detections above a S/N of 5.25. Be- 
low this S/N, the purity of the sample was found to 
decrease rapidly with only ~50% purity in the range 
5 < S/N < 5.25 (seven false detections). These results 
arc consistent with estimates of purity based on cross- 
identification of the ACT detections with other catalogs 
(See Section |4~1T ). From simple Gaussian statistics, one 
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expects fewer than five false detections in a sample se- 
lected with S/N > 5. Thus, some fraction of the false 
detections are may be the result of localized noise not 
accounted for by th e weighting and ma tched filter. 

Simulations from [Sehgal ct al.l (J2010D with noise from 
difference maps (see Section |2~3| were used to estimate 
completeness. Tabled] summarizes the findings. Due pri- 
marily to the uneven depth of coverage, the population 
of sources detected between 15 and 20 mJy was found 
on average to be 86% complete and the population be- 
tween 20 and 30 mJy to be 97% complete. Because of the 
strict 50 mJy lower bound set for detections in areas of 
the map characterized by integration times below 7 min- 
utes, the 50 - 90mJy range also suffers fro m an incom- 
pleten ess of 94%. The full simulations of iSchgal et al.l 
(|2010f ) include a correlation between radio sources and 
galaxy clusters, and we have excluded the SZ component 
in order to simplify the current study. At 148 GHz, clus- 
ters manifest themselves as arcminute-scale temperature 
decrements in the map which will cancel source flux in 
superposed source-cluster pairs. We also ra n the test de- 
scribe d here with the SZ component from iSehgal et al.l 
(|2010h included in the simulations. The resulting cancel- 
lation of source flux density by cluster decrements was 
found to cause an additional few percent of the incom- 
pleteness in the source population with flux densities be- 
low 30 mJy. 

4. DISCUSSION 

4.1. Comparison to Other Source Catalogs 

As a first step in ACT source characterization we con- 
sider cross-identifications with other catalogs. Matches 
are established within a 30" radius about an ACT source. 
The choice of association radius was made based on the 
positional rms of the ACT catalog and comparison cat- 
alogs (allowing for outliers) as well as the fact that the 
source of low-frequency radio signals in a given system 
may be physically displaced from the source of high- 
frequency radio signal. A general search through the 
NASA/IPAC Extragalactic Database finds that thirty- 
one ACT sources have measured redshifts, ranging from 
0.003 to 2.46. In what follows we consider in more detail 
several catalogs of particular relevance to the 148 GHz 
source population. 

Of our 157 sources, 109 match sources in the 5,8, and 
20 GHz AT20G catalog. There arc 180 AT20G sources in 
the survey area such that a random cross-identification 
would occur once in roughly 11600 cases H The AT20G 
catal og is incomplete below 100 mJy ([Murphy et al.l 
|2010[ ). Given that nearly all the radio sources detected 
at 1 48 GHz are expect ed to have relatively flat spectra 
fe.g.. lVieira et al. 2010), faint ACT sources may not have 
matches in AT20G (See Section l4~3l ). We have proposed 
for time on the Australia Telescope Compact Array to 
measure flux densities for sources in the ACT catalog 
that do not appear in AT20G. 

All but six of the ACT sources are co-located within 
30" of sources from AT20G or the 0.84 GHz Syd- 
ney University Mo longlo Sky Survey (SUMSS) catalog 
(|Mauch et al.ll200"3i) . Within our survey area, the sample 

4 This rough statistic estimates the probability of a spurious 
detection in the ACT data falling in the fractional area (TV X it x 
30" ) occupied by the N sources from the auxiliary catalog. 



of 14030 SUMSS sources is complete to 8 mJy. A random 
cross-identification with a SUMSS detection is a l-in-150 
event and thus a spurious SUMSS association is likely. 
Of the six ACT sources without cross-identification in 
AT20G or SUMSS, two (ACT-S J011830-511521, ACT- 
S J033133-515349) are within 50" of a SUMSS source, 
and the former is relatively bright with a 148 GHz flux 
density of 47.6 mJy. Furthermore, a preliminary reduc- 
tion of the ACT 218 GHz data identifies one of the re- 
maining four (ACT-S J031823-533148) as a ha detec- 
tion. ACT-S J034157-515140, ACT-S J004042-511830, 
and ACT-S J035343— 534553 have no match in the aux- 
iliary catalogs and may be false detections. This number 
of false detections is consistent with the study of sample 
purity presented in Section [3.31 

Comparing to the recently reported 2.0 mm measure- 
ment s from the South Pole Telescope (SPT) (jVieira et al.l 
I2010D . we find twenty- four c r oss-id entifications with ACT 
sources. The IVieira et al.l ([2010D study used a square 
survey of 87 square-degrees centered at 05h right ascen- 
sion. As such, the ACT and SPT surveys have only 
fractional overlap. Nevertheless, 2304 of the 3496 SPT 
source candidates (S/N > 3) fall within the ACT sur- 
ve y. All twenty-four matching sources were categorized 
in lVieira et al.l (|2010D as synchrotron-dominated. 

Finally, we compare the ACT catalog to the Infrared 
Astronomical Satellite Point Source Catalog (IRAS PSC; 
IHelou fc Walker] [l988t ). Three of the detections, ACT- 
S J041959-545622 (NGC 1566), ACT-S J04285-542959 
and ACT-S J033133-515352 (IC 1954), are identified 
with sources in the IRAS PSC. All three sources display 
lower frequency radio emission and have been identified 
in the preliminary ACT 218 GHz analysis. 

4.2. Correlation with X-ray Clusters 

Radio-loud AGN arc frequently found in Brightest 
Cluster Galaxies (BCGs). From a study of radio- loud 
AGNs in the SDSS using data from the National Ra- 
dio Astronomy Obs ervatory (NRAO) Ve ry Large Array 
(VLA) Sky Survey ([Condon et al.l 119981) and the Faint 
Imag es of the Radio Sk y At Twent y Cent imeters sur- 
vey ([Becker et al.l 119951) . iBest et al.l (|2007l ) found that 
the probability of a BCG hosting a radio-loud AGN is 
significantly enhanc ed compared to fiel d galaxies of the 
same stellar mass. iLin &: Mohrl (J2007I ) found that the 
fraction of BCGs being radio-loud is higher compared to 
that of cluster galaxies of similar luminosity. This en- 
hancement of radio activity is probably due to the fact 
that BCGs are located in special places - the centers of 
clusters - and the AGN activity is likely fueled by gas 
cooling, or due to galaxy interactions in these high den- 
sity regions. 

We performed a simple check to test whether any of 
the ACT 148 GHz sources were associated with clus- 
ters by cross matchin g against the REFLEX catalog 
(jBohringer et al.ll2004D . which is a homogeneous, X-ray 
selected sample with a nominal (0.1-2.4 keV) flux limit 
of > 3 x 10 -12 erg s _1 covering 5 < +2.5 dcg. The RE- 
FLEX catalog is > 90% complete above its nominal flux 
limit, and 23 REFLEX clusters are located within 2008 
ACT survey region, spanning the redshift range 0.03 < 
z < 0.34. Using a generous 5' matching radius, we find 
that three ACT sources are associated with REFLEX 
clusters: ACT-S J062142-524136 (RXC J0621. 7-5242); 
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Fig. 8. — Radio source spectral indices. In this radio color-color 
diagram, 5-20 GHz and 20-148 GHz spectral indices are shown for 
ACT-AT20G cross-identified sources. The population is dominated 
by sources which are peaked (lower right quadrant) or falling (lower 
left quadrant). Black (Gray) crosses correspond to sources with 
148 GHz flux density greater (less) than 50 mjy. The low flux 
sample is incomplete and suffers from selection bias that favors 
sources with more negative spectral indices. 



ACT-S J042906-534943 (RXC J0429.1-5350/AS0463); 
and ACT-S J062620-534136 (RXC J0626.3-5341/A3391). 
The separation between the ACT sources and the corre- 
sponding REFLEX cluster positions is 0.3 — 1.1'. All of 
these clusters are at very low redshift (0.041-0.055), and 
have low masses (~ (0.4 — 2) x 10 14 Mq, inferred from 
their X-ray luminosities). 

All three of these ACT sources have corresponding 
matches in the AT20G catalog. For ACT-S J062142- 
524136, both the ACT and AT20G sources are located 
within < 5" of each other and are coincident with the 
BCG, from inspection of DSS imaging. For the other 
two sources, we find that either the ACT source position 
(in the case of RXC J0626.3-5341) or the AT20G position 
(in the case of RXC J0429. 1-5350) is coincident with the 
BCG, although the AT20G and ACT positions are off- 
set by < 19". In all cases, the projected radial distance 
between the BCG and the REFLEX X-ray position is 
< 60 kpc. 

4.3. Source Spectra 

It is established that radio source spectra extending to 
148 GHz are not well characteriz e d by a si mple power law 
S(v) e x v a (e.g., Ide Zotti et ail ([2010l )). iMurphv et all 
(|2010T l used a color-color comparison of spectral indices 
at 5-8 GHz and 8-20 GHz to show that the AT20G pop- 
ulation may be decomposed into classical steep (and 
steepening) spectrum sources, sources that peak between 
their bands, and sources that show fiat, rising or up- 
turned spectra. Following this example we construct 
a color-color comparison of 5-20 GHz and 20-148 GHz 
measurements, where the 5 and 20 GHz flux densities 
are from AT20G. The variability of these sources makes 
a per-source comparison difficult, but a study using the 
109 ACT-AT20G cross-identifications as an ensemble is 
meaningful. Figure [8] is the as_2o vs. «2o-i48 color- 
color diagram. The figure shows that the ACT-AT20G 



cross-identified sources are predominantly characterized 
by spectral steepening. The sample is biased towards 
steepened spectra at low flux density (Si48 < 50 mjy) 
due to incompleteness in AT20G. 

The population can be further divided according to 
148 GHz flux density. In Figure the black crosses cor- 
respond to the forty-two brightest sources in the ACT 
sample, and the gray crosses correspond to the faint half. 
The dividing flux density, 50 mjy, was chosen such that 
all but two of the ACT detections in the brighter sample 
have cross-identifications in AT20G. As described in Sec- 
tion al) below this flux density the mean spectral indices 
of the population of ACT- AT20G cross- identified sources 
is biased negative by the incompleteness in AT20G be- 
low 100 mjy. The two sources in this high-flux subset 
which do not have AT20G counterparts are character- 
ized by 150 GHz flux densities close to the 50 mjy cut- 
off and may have flatter spectra than the sample aver- 
age. Including two extra sources with a = changes 
the mean spectral indices by ~ 25% of the statistical er- 
ror, and therefore we simply use the forty sources with 
AT20G counterparts in the following analysis. Consid- 
ering only the unbiased, bright half of the distribution, 
the average spectrum steepens between 5-20 GHz and 
20-148 GHz. The median spectral indices of the unbi- 
ased sample of sources are 05-20 = —0.07 ± 0.37(±0.06) 
and a 20 -i48 = -0.39 ± O.24(±O.O4)0 In obtaining the 
spectral indices we compare the deboosted flux densi- 
ties from ACT with the raw flux densities from AT20G. 
Using the raw AT20G should not significantly bias the 
index estimates because the AT20G detections are all 
characterized by a S/N greater than 15. 



0.2 



T 1 




+ + + 



10" 



148 GHz Flux fJvl 



Fig. 9. — Radio color-magnitude diagram. As in Figure [8] 
the data are divided between high flux densities (black crosses) 
and low flux densities (gray crosses) at 50 mjy. The unbiased 
50 mjy sample has a 5-148 GHz spectral index of 05-148 = 
-0.20 ± 0.21 (±0.03). The low flux density data suffer a selection 
bias that excludes sources with flat or rising spectra. 



When restricted to the unbiased sample with flux den- 
sities above 50 mjy, the 5-148 GHz spectral index is 
a 5 _i 48 = -0.20 ± 0.21(±0.03). This distribution is in 

5 The errors on spectral indices are the 68% confidence levels of 
the distribution and, in parentheses, for the median. 
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2.5 a tension with the SPT-reported mean 5-150 GHz 
spectral slop e (for 57 sources) of 0:5-150 = —0.13 ± 
0.21(±0.03) (jVieira et al.l l2010h FI IVieira et all (poTTih 
claim that the mean spectral index of the synchrotron- 
dominated species remains near —0.1 to 2.0 mm (w 
148 GHz) after which it steepens such that the aver- 
age slope between 2 mm and 1.4 mm (« 220 GHz) is 
—0.5. This study suggests that the transition to the steep 
spectrum is more gradual. In fact, the spectral slope 
«2o-i48 = — 0.39±0.24(±0.04) approaches t he —0.5 slope 
betwe en 2.0 mm and 1.4 mm reported in IVieira et al.l 
(|2010[ ). This picture is further supported by the rescaling 
of ACT 148 GHz source counts relative to source counts 
at 30 GHz (See Section EU). Figure M shows the ACT 
5-148 GHz spectral indices as a function of flux density. 
The low flux density sample, represented by gray crosses, 
is incomplete for high spectral indices. Follow-up of the 
148 GHz selected sources without matches in AT20G will 
complete the picture in the range 20-50 mJy. 

4.4. Source Counts 

The differential number counts for ACT sources based 
on data in Table Q] are plotted in Figure [TUJ The figure 
shows that the ACT counts are fit reasonably well by 
the model for radio sources from iToffolatti et al.l (|1998[ ) 
scaled by a factor of 0.34 ± 0.04. Depending on whether 
the rescaling is fit to the data in log or linear coordinates 
the best fit value varies from 0.31 to 0.35. As in Fig- 
ure [TU1 we adopt the rescaling of 0.34 which lies between 
these two values. The ACT res ults are also consi stent 
with counts reported by the SPT (IVieira et al.H2010D. As 
part of the WMAP 7-year analysis. I Wright et all (J2009D 
fit di fferential source cou nts at 3 GHz from WMA P, the 
VSA (IClearv et al.lj200fih CBI UMason et all 120031) . and 
DASI (iKovac et all 120021) . They fou nd that the best fit 
scaling for the IToffolatti et al.l (J1998I ) model at 30 GHz is 
0.64. This result was also consistent with 30 GHz mea- 
surements by the SZA which covere d a fainter range of 
flux densities (|Muchovei et alll2010D . The apparent dis- 
crepancy between the 30 GHz and 148 GHz scalings can 
be explained by the steepe ning of the radio source spec- 
trum described in Section J4.3I The source coun t model 
traditionally adopted from lToffolatti et all (|1998l) uses an 
average spectral index of a = for radio sources between 
20 and 200 GHz. From the unbiased sample of spectral 
indices (dark crosses) in Figure the slope between 20 
and 148 GHz is a re —0.4. As can be seen from Figure [TUl 
the radio source number counts are well-approximated by 
a power law N(>S) oc S _1 between 0.01 and 1 Jy (dN/dS 
S 5/2 K gi/2)_ It f n ows t h a t the 148 GHz number counts 
should be rescaled by 



0.64 x 



5*148 



GHz 



'30 GHz 



/l48GHz\ Q 
! °- 64X UgHz-) (9) 

;0.64x 5~ a4 
;0.35, 



consistent with the scaling of IToffolatti et ail (J1998I) to 
the ACT data. Also shown in the plot are radio source 

6 Note the tension arises for the center of the spectral index 
distribution which is better constrained than the index of any given 
source. 



models from |Pe Zotti et all (120051) and iSehgal et all 
(|2010l) . The iDe Zotti et al.l (|2005l ) model is consistent 
with the ACT source counts at low flux densities (S < 0.1 
Jy) and over-predicts the counts at higher flu x densities. 
The model presented in ISehgal et all (J2010I ) , although 
underestimating the source counts at S < 0.1 Jy, seems 
to be consistent with ACT data at higher fluxes. 

Differential source c ounts derived from t he 150 GHz 
catalog fr om ACBAR flR cichardt et al. 2009) are consis- 
tent with IToffolatti et all (|1998l ) scaled by 0.64. How- 
ever, ACBAR's resolution and sensitivity made it sen- 
sitive to only the brightest sources with flux densities 
mainly in excess of 100 mJy where there are fewer 
sources. Furthermore, the ACBAR fields were chosen 
specifically to ha ve bright quasars for beam measurement 
(|Kuo et al.ll2007l ) such that the differential source counts 
derived from the ACBAR data are biased high. The ACT 
field was chosen with no a-priori knowledge of source 
population. Allowing for variability, there is reasonable 
agreement in the flux densities reported by ACBAR and 
ACT for the seven sources common to both catalogs. 

Also shown in Figure [TOl ar c models for counts of dusty 
starb u rst galaxies i Toffolatti et al. 1998; Lag ache et all 
[200l iNegrelloet all 120071 ). The brightest infrared 
sources in these models are 10 mJy. Given that all 
sources in the ACT catalog have flux densities greater 
than 10 mJy, these models predict that the 148 GHz 
selected ACT catalog should have few or no infrared 



4.5. Contribution to the Power Spectrum 

The ^=3000-10000 power spectrum of the 148 GHz 
sky is dominated by synchrotron and infrared sources 
as well as the ther mal SZ from clusters (jHall et alll2010l : 
iFowler et al.ll2010j ) . In particular, the power spectrum at 
the highest multipoles constrains the Poisson distributed 
component of the source population, and an understand- 
ing of the residual synchrotron population helps break 
the high-^ spectral degeneracy between synchrotron and 
infrared sources. The contribution to the power spec- 
trum by a Poisson-distributed population of sources is a 
function of the number counts: 



C 



PS 



6B„ 

ST 



"-*>■ 



(10) 



where Su m is the upper-limiting flux density of the resid- 
ual (i.e., unmasked) sources in the data. With a lim- 
iting flux densi t y of 2 mJy and the rescaled model of 
IToffolatti et al.l (|1998l) from Section 14.41 one expects a 
synchrotron contribution to the Poisson power spectrum 
o f C Sync = (2.8 ±0.3) x 10" VK 2 . 

IFowler et all (J2010D used a 20 mJy cut for masking 
sources and found a Poisson spectrum from all residual 
sources of C PS = (7.8 ± 2.3) x 10~VK 2 . Thus from this 
study we expect an infrared contribution to the Poisson 
spectrum of C IR = (5.0 ± 2.3) x lO^uK 2 . A s imilar 
argument applied to the study of iLueker et al.l (|2010f ) 
results in an estimate of the residual Poisson term for 
infrared so urces of C m » (6.3 ± 0.5) x 10~VK 2 . Thus, 
within the Fowler e t al.l (J2010D errors, the two studies of 
the high-^ power spectrum at 148 GHz are consistent. 

5. CONCLUSIONS 
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Fig. 10. — 148 GHz Differential Source Counts. Derived from Table[JJand corrected for completeness, the ACT differential source counts 
(bold diamonds) arc plotted together with models of radio and infrared source populations at 143 GHz (as computed for the 2 mm band of 
the P lanck satellite). The data are consistent with being dominated by radio sources. Data (gray diamonds) from the SPT l|Vieira et al.l 
2010) arc consistent with the ACT counts. Errors (lcr) are Poissonian. 



We have presented results on extragalactic sources at 
148 GHz from data taken by ACT during the 2008 ob- 
serving season. A catalog of 157 millimeter sources has 
been presented with sources detected across two decades 
in flux density, from 15mJy to 1500 mjy. The flux den- 
sity calibration of the sources derives from observations 
of Uranus with 6% error. Bias in the quoted flux densi- 
ties due to beam shape uncertainty is estimated at less 
than 1%. Typical statistical lcr errors for the source flux 
density range from 2.5 mjy to 5 mjy. The catalog as- 
trometry error for the brightest sources is characterized 
by an rms of 3.5". The catalog is estimated to be 98% 
pure and complete above 20 mjy. 

Comparison to other catalogs shows that 98% of the 
ACT detections correspond to sources detected at lower 
radio frequencies. The differential source counts are also 
consistent with the finding that ACT detections corre- 
spond to sources detected at lower radio frequencies. 
In particular, the sour ce counts are fi t reaso nably well 
by the radio model of Toffolatt i et al.l (I1998D scaled by 
0.34. This scaling, compared t o a scaling 0.64 found at 
30 GHz bv IWright et all (|2009f ). suggests that the pop- 
ulation of radio sources is characterized, on average, by 
spectral steepening between 30 GHz and 148 GHz. This 
conclusion is consistent with the average spectral indices 
derived from the combined AT20G and ACT datasets. 
Future work will addres s the m ore involved comparison 
with the iDe Zotti et all (|2005h and iSehgal et all (|2010h 
radio source models. With the rescaled model from 
iToffolatti et al.l (1 19981 ) and a 20 mjy cut, the residual 
contribution of the synchrotron population to the Pois- 
son power spectrum is C Sync = (2.8 ± 0.3) x 10" 6 u,K 2 . 



Future ACT source work will incorporate the 218 GHz 
and 277 GHz bands, deeper coverage integrating the 
2007, 2009, and (ongoing) 2010 seasons, as well as the 
equatorial survey overlapping the deep Sloan Digital Sky 
Survey Stripe 82. 



The ACT project was proposed in 2000 and funded 
by the U.S. National Science Foundation on January 1, 
2004. Many have contributed to the project since its in- 
ception. We especially wish to thank Asad Aboobaker, 
Christine Allen, Dominic Benford, Paul Bode, Kristen 
Burgess, Angelica de Oliveira-Costa, Peter Hargrave, 
Norm Jarosik, Amber Miller, Carl Reintsema, Felipe Ro- 
jas, Uros Seljak, Martin Spergel, Johannes Staghun, Carl 
Stahle, Max Tegmark, Masao Uehara, Katerina Visnjic, 
and Ed Wishnow. It is a pleasure to acknowledge Bob 
Margolis, ACT's project manager. Reed Plimpton and 
David Jacobson worked at the telescope during the 2008 
season. ACT is on the Chajnantor Science preserve, 
which was made possible by the Chilean Comision Na- 
tional de Investigation Cientifica y Tecnologica. 

This work was supported by the U.S. National Science 
Foundation through awards AST-0408698 for the ACT 
project, and PHY-0355328, AST-0707731 and PIRE- 
0507768. Funding was also provided by Princeton Uni- 
versity and the University of Pennsylvania. The PIRE 
program made possible exchanges between Chile, South 
Africa, Spain and the US that enabled this research pro- 
gram. Computations were performed on the GPC su- 
percomputer at the SciNet HPC Consortium. SciNet is 
funded by: the Canada Foundation for Innovation under 



14 



the auspices of Compute Canada; the Government of On- 
tario; Ontario Research Fund - Research Excellence; and 
the University of Toronto. 

TM was supported through NASA grant 
NNX08AH30G. JBJ was supported by the FONDECYT 
grant 3085031. ADH received additional support from 
a Natural Science and Engineering Research Council 
of Canada (NSERC) PGS-D scholarship. AK and BP 
were partially supported through NSF AST-0546035 
and AST-0606975, respectively, for work on ACT. HQ 
and LI acknowledge partial support from FONDAP 



Centro de Astrofisica. NS is supported by the U.S. 
Department of Energy contract to SLAC no. DE- 
AC3-76SF00515. RD was supported by CONICYT, 
MECESUP, and Fundacion Andes. RH was supported 
by the Rhodes Trust. ES acknowledges support by 
NSF Physics Frontier Center grant PHY-01 14422 to 
the Kavli Institute of Cosmological Physics. YTL 
acknowledges support from the World Premier Inter- 
national Research Center Initiative, MEXT, Japan. 
The ACT data will be made public through LAMBDA 
(http:/ /lambda. gsf c.nasa.gov/[) and the ACT website 
( http://www.physics.princeton.edu/act/[ ). 



APPENDIX 



REFERENCES 



Austcrmann, J. E., et al. 2010, MNRAS, 401, 160 
Battaglia, N., Bond, J. R., Pfrommer, C, Sievers, J. L., & 

Sijacki, D. 2010, ApJ, 725, 91 
Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 
Best, P. N., von der Linden, A., Kauffmann, G., Heckman, T. M., 

& Kaiser, C. R. 2007, MNRAS, 379, 894 
Bohringer, H., et al. 2004, A&A, 425, 367 

Cleary, K. A., et al. 2006, VizieR Online Data Catalog, 736, 340 
Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., Perley, 

R. A., Taylor, G. B., & Broderick, J. J. 1998, AJ, 115, 1693 
Coppin, K., Halpern, M., Scott, D., Borys, C, & Chapman, S. 

2005, MNRAS, 357, 1022 
Crawford, T. M., Switzer, E. R., Holzapfel, W. L., Reichardt, 

C. L., Marrone, D. P., & Vieira, J. D. 2010, ApJ, 718, 513 
de Zotti, G., Massardi, M., Negrello, M., & Wall, J. 2010, 

A&A Rev., 18, 1 
De Zotti, G., Ricci, R., Mesa, D., Silva, L., Mazzotta, P., 

Toffolatti, L., & Gonzalez-Nuevo, J. 2005, A&A, 431, 893 
Fowler, J. W., et al. 2010, ApJ, 722, 1148 
Hall, N. R., et al. 2010, ApJ, 718, 632 
Helou, G. & Walker, D. W., eds. 1988, Infrared astronomical 

satellite (IRAS) catalogs and atlases. Volume 7: The small scale 

structure catalog, Vol. 7 
Hincks, A. D., et al. 2009, ApJS, 191, 423 
Hogbom, J. A. 1974, A&AS, 15, 417 
Kesteven, M. J. L., Bridle, A. H., & Brandie, G. W. 1977, AJ, 82, 

541 
Kovac, J. M., Leitch, E. M., Prykc, C, Carlstrom, J. E., 

Halverson, N. W., & Holzapfel, W. L. 2002, Nature, 420, 772 
Kuo, C. L., et al. 2007, ApJ, 664, 687 
Lagache, G., et al. 2004, ApJS, 154, 112 

Lima, M., Jain, B., Devlin, M., & Aguirre, J. 2010, ApJ, 717, L31 
Lin, Y. & Mohr, J. J. 2007, ApJS, 170, 71 
Lin, Y., Partridge, B., Pober, J. C, Bouchefry, K. E., Burke, S., 

Klein, J. N., Coish, J. W., & Huffenberger, K. M. 2009, ApJ, 

694, 992 
Luckcr, M., et al. 2010, ApJ, 719, 1045 



Mason, B. S., et al. 2003, ApJ, 591, 540 

Mauch, T., Murphy, T., Buttery, H. J., Curran, J., Hunstead, 

R. W., Piestrzynski, B., Robertson, J. C, & Sadler, E. M. 

2003, MNRAS, 342, 1117 
Melin, J., Bartlett, J. G., & Delabrouille, J. 2006, A&A, 459, 341 
Menanteau, F., et al. 2010, ApJ, 723, 1523 
Muchovej, S., et al. 2010, ApJ, 716, 521 
Murphy, T., et al. 2010, MNRAS, 402, 2403 
Negrello, M., Perrotta, F., Gonzalez-Nuevo, J., Silva, L., De 

Zotti, G., Granato, G. L., Baccigalupi, C, & Danese, L. 2007, 

MNRAS, 377, 1557 
Nolta, M. R., et al. 2009, ApJS, 180, 296 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, 

B. P. 1992, Numerical Recipes in C: The Art of Scientific 

Computing, 2nd cdn. (Cambridge University Press) 
Reichardt, C. L., et al. 2009, ApJ, 701, 1958 
Scott, K. S., et al. 2008, MNRAS, 385, 2225 
Sehgal, N., et al. 2010, ApJ, 709, 920 
Shaw, L. D., Nagai, D., Bhattacharya, S., & Lau, E. T. 2010, 

ApJ, 725, 1452 
Staniszewski, Z., et al. 2009, ApJ, 701, 32 
Sunyaev, R. A. & Zel'dovich, Y. B. 1970, Comments on 

Astrophysics and Space Physics, 2, 66 
Swetz, D. S., et al. 2010, arXiv: 1007.0290 
Tegmark, M. & de Oliveira-Costa, A. 1998, ApJ, 500, L83 
Toffolatti, L., Argueso Gomez, F., De Zotti, G., Mazzei, P., 

Franceschini, A., Danese, L., & Burigana, C. 1998, MNRAS, 

297, 117 
Trac, H., Bode, P., & Ostriker, J. P. 2011, ApJ, 727, 94 
Valtaoja, E., Terasranta, H., Urpo, S., Nesterov, N. S., Lainela, 

M., & Valtonen, M. 1992, A&A, 254, 80 
Vanderlinde, K., et al. 2010, ApJ, 722, 1180 
Vieira, J. D., et al. 2010, ApJ, 719, 763 
Weiland, J. L., et al. 2011, ApJS, 192, 19 
Wright, E. L., et al. 2009, ApJS, 180, 283 



TABLE Al 
ACT High Significance 148 GHz Extragalactic Source Catalog 
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ACT ID 



RA (J2000) Dec 
h m s ° ' 



S/N 



o a 

(mjy) 



t>db 

(mjy) 



tint C 

(min) 



AT20G ID 



ACT-S J001850-511454 


00 


18 


50 


-51 


14 


51 


10.4 


53.7 


53.7±|;1 


ACT-S J001851-492949 


00 


18 


51 


-49 


29 


19 


7.1 


52.8 


52.8j£j 


ACT-S J002513-542737 


00 


25 


13 


-54 


27 


37 


6.7 


26.1 


2A.8ti 9 


ACT-S J003134-514308 


00 


31 


34 


-51 


43 


08 


13.9 


52.7 


r:o 7+3.8 
oz, '-3.8 


ACT-S J003734-530729 


00 


37 


34 


-53 


07 


29 


11.9 


39.1 


QQ r + 3.3 

38.t>_3.3 


ACT-S J004012-514950 


00 


40 


12 


-51 


49 


50 


5.7 


20.2 


18.71^ 


ACT-S J004042-511830 


00 


40 


12 


-51 


18 


30 


5.6 


22.9 


21 1+ 4 - 2 


ACT-S J004906-552106 


00 


49 


06 


-55 


21 


06 


14.1 


62.6 


62.6l*j 


ACT-S J004949-540240 


00 


49 


19 


-54 


02 


10 


6.7 


21.8 


20.7±1J 


ACT-S J005240-531127 


00 


52 


40 


-53 


11 


27 


5.5 


17.6 


16 2+ 3 ' 4 


ACT-S J005605-524154 


00 


56 


05 


-52 


41 


51 


8.2 


28.6 


27.6±H 


ACT-S J005622-531845 


00 


56 


22 


-53 


18 


15 


9.7 


32.4 


3i-6±|:! 


ACT-S J005706-521423 


00 


57 


06 


-52 


14 


23 


13.9 


47.4 


46.8i 3 3 ;l 


ACT-S J005855-521925 


00 


58 


55 


-52 


19 


25 


20.3 


70.4 


70.4T.|i 


ACT-S J010309-510907 


01 


03 


09 


-51 


09 


07 


8.8 


33.5 


qo c+3-8 
oi - a -3.8 


ACT-S J010330-513544 


01 


03 


30 


-51 


35 


11 


11.7 


41.9 


41 0+3-6 


ACT-S J011324-532938 


01 


13 


21 


-53 


29 


38 


7.2 


22.4 


2i-4±|;l 


ACT-S J011654-544653 


01 


16 


51 


-54 


46 


53 


5.8 


21.0 


19.511;* 


ACT-S J011830-511521 


01 


18 


30 


-51 


15 


21 


13.0 


48.4 


47.7+lf 


ACT-S J011950-535714 


01 


19 


50 


-53 


57 


11 


46.9 


157.2 


157.2li 4 


ACT-S J012006-521104 


01 


20 


06 


-52 


11 


01 


6.1 


19.6 


18.4±!;! 


ACT-S J012457-511309 


01 


24 


57 


-51 


13 


09 


51.0 


199.1 


I99.l±l;g 


ACT-S J012501-532514 


01 


25 


01 


-53 


25 


11 


5.3 


15.5 


14.2i 3 3 ;i. 


ACT-S J012623-510305 


01 


26 


23 


-51 


03 


05 


7.0 


26.0 


24.811:1 


ACT-S J012755-513642 


01 


27 


55 


-51 


36 


12 


12.1 


39.5 


qo q+3.3 
38.8_ 33 


ACT-S J012756-532933 


01 


27 


56 


-53 


29 


33 


9.6 


27.6 


26.9+f° 


ACT-S J012834-525520 


01 


28 


31 


-52 


55 


20 


6.3 


18.5 


174+3:0 


ACT-S J013107-545707 


01 


31 


07 


-54 


57 


07 


6.3 


23.8 


00 q+3.9 
zz '°-3.9 


ACT-S J013225-512903 


01 


32 


25 


-51 


29 


03 


5.3 


18.4 


16.811;? 


ACT-S J013306-520003 


01 


33 


06 


-52 


00 


03 


74.7 


239.0 


239.0+!'2 


ACT-S J013409-552612 


01 


34 


09 


-55 


26 


12 


6.0 


24.7 


23 l+ 4 -' 2 


ACT-S J013541-514943 


01 


35 


11 


-51 


49 


13 


6.6 


19.8 


18.8Ip 


ACT-S J013548-524414 


01 


35 


18 


-52 


44 


11 


18.2 


51.6 


Bi-6±l:l 


ACT-S J013727-543942 


01 


37 


27 


-54 


39 


12 


7.1 


25.3 


24-2±|;| 


ACT-S J013949-521739 


01 


39 


19 


-52 


17 


39 


12.1 


35.2 


34.6±i| 


ACT-S J014648-520232 


01 


46 


18 


-52 


02 


32 


25.3 


77.5 


77.5ll;i 


ACT-S J015358-540649 


01 


53 


58 


-54 


06 


49 


11.4 


33.7 


QQ 1+3.0 

• 3,3 - L -3.0 


ACT-S J015420-510750 


01 


54 


20 


-51 


07 


50 


49.3 


172.1 


m.irj-j 


ACT-S J015559-512538 


01 


55 


59 


-51 


25 


38 


5.8 


18.2 


16.91!;! 


ACT-S J015649-543940 


01 


56 


19 


-54 


39 


40 


13.4 


42.2 


41.6T.1J 


ACT-S J015817-500415 


01 


58 


17 


-50 


04 


15 


5.8 


24.6 


22.8T 4 ; 4 


ACT-S J015914-530902 


01 


59 


11 


-53 


09 


02 


9.2 


25.3 


24.712-8 


ACT-S J020448-550257 


02 


04 


18 


-55 


02 


57 


14.7 


48.4 


47.911;! 


ACT-S J020649-534528 


02 


06 


19 


-53 


45 


28 


5.1 


15.2 


13.91!;° 


ACT-S J020920-522921 


02 


09 


20 


-52 


29 


21 


11.1 


30.5 


29.9ll:| 


ACT-S J021046-510100 


02 


10 


46 


-51 


01 


00 


470.3 


1677.6 


1678.0l!'g 


ACT-S J021519-510435 


02 


15 


19 


-51 


04 


35 


5.5 


17.3 


15.9l!i 


ACT-S J021603-520007 


02 


16 


03 


-52 


00 


07 


12.7 


34.6 


34.il 2 2 ; 7 r 


ACT-S J021709-542750 


02 


17 


09 


-54 


27 


50 


5.8 


17.3 


16.llli 


ACT-S J021835-550354 


02 


18 


35 


-55 


03 


51 


12.5 


43.0 


42.4l!;l 


ACT-S J022216-510627 


02 


22 


16 


-51 


06 


27 


9.1 


29.7 


28.91!;^ 


ACT-S J022330-534737 


02 


23 


30 


-53 


47 


37 


40.5 


119.9 


119.91!;° 


ACT-S J022530-522547 


02 


25 


30 


-52 


25 


17 


15.3 


40.2 


39.8l2;f 


ACT-S J022820-553732 


02 


28 


20 


-55 


37 


32 


5.6 


21.9 


20.2 +4: J 


ACT-S J022821-554601 


02 


28 


21 


-55 


46 


01 


25.9 


104.7 


104.711°, 


ACT-S J022912-540324 


02 


29 


12 


-54 


03 


21 


112.5 


333.1 


333.111;° 


ACT-S J022925-523226 


02 


29 


25 


-52 


32 


26 


17.2 


45.3 


44.9l 2 2 j 



6.1 
3.0 
10.0 
10.3 
13.9 
12.2 
9.5 
7.5 
14.2 
14.5 
12.5 
13.5 
12.7 
12.4 
10.4 
12.0 
15.4 
11.3 
11.0 
15.4 
14.8 
11.4 
17.6 
11.1 
14.1 
18.1 
17.7 
10.3 
12.9 
17.1 
8.8 
16.7 
18.8 
12.6 
17.7 
16.0 
16.9 
14.2 
15.0 
15.1 
8.3 
19.8 
13.6 
19.4 
20.2 
13.6 
15.1 
20.1 
16.4 
12.5 
15.1 
19.4 
21.4 
9.8 
9.0 
19.9 
21.7 



ATG20 J001849-511455 

ATG20 J002511-542739 
ATG20 J003134-514325 
ATG20 J003735-530733 



ATG20 J004905-552110 



ATG20 J005705-521418 
ATG20 J005855-521926 
ATG20 J010306-510907 
ATG20 J010329-513551 
ATG20 J011323-532949 



ATG20 J011950-535717 
ATG20 J012008-521102 
ATG20 J012457-511316 

ATG20 J012624-510309 
ATG20 J012756-513641 

ATG20 J012834-525519 



ATG20 J013305-520003 
ATG20 J013408-552616 
ATG20 J013540-514945 
ATG20 J013548-524417 



ATG20 
ATG20 
ATG20 
ATG20 
ATG20 
ATG20 
ATG20 
ATG20 



J013949- 
J014648- 
J015358- 
J015419- 
J015557- 
J015649- 
J015817- 
J015913- 



521746 
520233 
540653 
510751 
512545 
543949 
500419 
530853 



ATG20 J020647-534543 



ATG20 J021046-510101 



ATG20 J021603-520012 



ATG20 
ATG20 
ATG20 
ATG20 
ATG20 
ATG20 
ATG20 
ATG20 



J021834- 
J022215- 
J022330- 
J022529- 
J022820- 
J022821- 
J022912- 
J022925- 



550350 
510629 
534740 
522555 
553725 
554603 
540324 
523226 
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RA (J2000) Dec 



S/N 



s a 

(mjy) 



Sdb 

(mjy) 



tint C 

(min) 



AT20G ID 



ACT-S J023245-535634 


02 


32 


45 


-53 


56 


31 


8.1 


21.9 


21.2+. 2 ; 7 . 


ACT-S J023357-503014 


02 


33 


57 


-50 


30 


11 


6.3 


21.4 


20.2+H 


ACT-S J023924-510824 


02 


39 


21 


-51 


08 


21 


8.5 


25.2 


24 4+3.0 
zqi -^-3.0 


ACT-S J024040-542942 


02 


40 


40 


-54 


29 


12 


6.0 


21.7 


20.7JtJ 


ACT-S J024154-534541 


02 


41 


51 


-53 


45 


11 


6.7 


17.8 


ie.9tll 


ACT-S J024313-510512 


02 


43 


13 


-51 


05 


12 


11.5 


35.1 


34.5i!;l 


ACT-S J024539-525756 


02 


45 


39 


-52 


57 


56 


7.6 


19.1 


18.4+1J 


ACT-S J024614-495346 


02 


46 


11 


-49 


53 


46 


10.3 


42.8 


41.9TJ;! 


ACT-S J024947-555619 


02 


49 


47 


-55 


56 


19 


5.1 


22.8 


20.9JH 


ACT-S J025112-520822 


02 


51 


12 


-52 


08 


22 


5.7 


15.0 


13 9+ 2 ' 7 


ACT-S J025204-514554 


02 


52 


01 


-51 


45 


51 


5.4 


14.4 


i n 9 + 2.8 
±0 - -2.8 


ACT-S J025328-544151 


02 


53 


28 


-54 


41 


51 


246.7 


837.5 


837.5+1* 


ACT-S J025629-522711 


02 


56 


20 


-52 


27 


11 


5.3 


14.1 


12.8+2'g 


ACT-S J025838-505200 


02 


58 


38 


-50 


52 


00 


33.2 


113.6 


113.6±|-* 


ACT-S J025849-533210 


02 


58 


19 


-53 


32 


10 


7.2 


19.5 


18.6j£j 


ACT-S J030056-510229 


03 


00 


56 


-51 


02 


29 


11.2 


36.8 


36.1+jj 


ACT-S J030132-525603 


03 


01 


32 


-52 


56 


03 


15.0 


39.7 


39 s+ 2 ' 6 


ACT-S J030327-523427 


03 


03 


27 


-52 


34 


27 


24.8 


66.4 


66.4+1 7 


ACT-S J030616-552808 


03 


06 


16 


-55 


28 


08 


15.2 


56.4 


56.5+3-7 


ACT-S J031206-554135 


03 


12 


06 


-55 


41 


35 


18.3 


72.7 


79 7+4:0 

,z -'-4.0 


ACT-S J031426-510430 


03 


14 


26 


-51 


04 


30 


29.0 


90.6 


90-5+11 


ACT-S J031823-533148 


03 


18 


23 


-53 


31 


18 


18.6 


51.6 


51.6+1* 


ACT-S J031910-500031 


03 


19 


10 


-50 


00 


31 


7.5 


30.2 


28.9±*-J 


ACT-S J032207-535419 


03 


22 


07 


-53 


54 


19 


6.5 


18.3 


-,70+2:9 


ACT-S J032212-504229 


03 


22 


12 


-50 


42 


29 


11.0 


35.5 


34.8+1? 


ACT-S J032327-522627 


03 


23 


27 


-52 


26 


27 


14.9 


41.0 


40.5+1* 


ACT-S J032650-533658 


03 


26 


50 


-53 


36 


58 


17.0 


48.0 


47.6+.2J 


ACT-S J033002-503518 


03 


30 


02 


-50 


35 


18 


6.5 


21.9 


20.7+j-j 


ACT-S J033114-524149 


03 


31 


11 


-52 


41 


49 


9.0 


24.4 


23.7+1? 


ACT-S J033126-525829 


03 


31 


26 


-52 


58 


29 


11.9 


31.6 


31.0+1 ?, 


ACT-S J033133-515352 


03 


31 


33 


-51 


53 


52 


5.9 


16.1 


15.0+1 8 , 


ACT-S J033444-521851 


03 


34 


11 


-52 


18 


51 


7.1 


19.2 


18.4i|l 


ACT-S J033554-543028 


03 


35 


51 


-54 


30 


28 


14.8 


44.0 


43.6+1° 


ACT-S J034157-515140 


03 


41 


57 


-51 


51 


10 


6.3 


16.6 


15.7+1 7 7 


ACT-S J034348-524112 


03 


43 


18 


-52 


41 


12 


8.3 


21.9 


21 2+ 2 ' 7 


ACT-S J034940-540111 


03 


49 


10 


-54 


01 


11 


7.1 


21.7 


20.7+1°, 


ACT-S J035128-514254 


03 


51 


28 


-51 


42 


51 


50.7 


145.2 


145.2+ 2 -9 


ACT-S J035343-534553 


03 


53 


13 


-53 


45 


53 


5.6 


16.3 


15.0+1° 


ACT-S J035700-495549 


03 


57 


00 


-49 


55 


19 


14.1 


58.3 


58.3_ 4 ]j 


ACT-S J035840-543403 


03 


58 


10 


-54 


34 


03 


8.3 


27.4 


26.4+3 1 '! 


ACT-S J040401-552022 


04 


04 


01 


-55 


20 


22 


7.0 


27.6 


26.3+1°, 


ACT-S J040621-503509 


01 


06 


21 


-50 


35 


09 


5.5 


20.5 


18.8+lt 


ACT-S J041137-514918 


04 


11 


37 


-51 


49 


18 


22.8 


72.5 


72.5+1?, 


ACT-S J041247-560044 


04 


12 


17 


-56 


00 


11 


5.1 


23.4 


21-4+1J 


ACT-S J041313-533157 


04 


13 


13 


-53 


31 


57 


13.2 


39.5 


38-9+1°, 


ACT-S J041959-545618 


04 


19 


59 


-54 


56 


18 


7.2 


23.8 


99 7+3.4 

zz -'-3.4 


ACT-S J042503-533201 


04 


25 


03 


-53 


32 


01 


49.8 


152.1 


152.1+lJ 


ACT-S J042842-500532 


04 


28 


12 


-50 


05 


32 


33.1 


154.7 


154.7+1 7 , 


ACT-S J042852-543001 


04 


28 


52 


-54 


30 


01 


6.8 


21.8 


20-8+13 


ACT-S J042906-534943 


04 


29 


06 


-53 


49 


13 


31.7 


87.4 


87.4+1! 


ACT-S J043221-510926 


04 


32 


21 


-51 


09 


26 


26.5 


86.8 


86.8+13 


ACT-S J043651 -521632 


04 


36 


51 


-52 


16 


32 


14.1 


38.4 


38.0+1 7 . 


ACT-S J044115-543848 


04 


41 


15 


-54 


38 


18 


6.1 


19.4 


18.2+13 


ACT-S J044158-515456 


04 


41 


58 


-51 


54 


56 


30.4 


87.5 


87.5+JI 


ACT-S J044502-523426 


04 


45 


02 


-52 


34 


26 


9.1 


25.5 


24-8+1 8 , 


ACT-S J044746-515103 


04 


47 


16 


-51 


51 


03 


9.7 


27.2 


26.5+1! 


ACT-S J044821-504140 


04 


48 


21 


-50 


41 


40 


10.8 


37.1 


36-3+14 


ACT-S J045028-534659 


04 


50 


28 


-53 


46 


59 


8.9 


24.8 


94 1+2-8 



20.7 
13.0 
16.8 
15.3 
20.9 
16.2 
23.4 
8.5 
8.2 
21.8 
20.7 
15.3 
21.0 
14.6 
20.6 
14.0 
21.6 
21.0 
10.9 
9.5 
15.2 
19.6 
9.3 
19.6 
14.5 
19.7 
18.9 
13.5 
21.1 
21.5 
19.8 
20.4 
16.6 
22.1 
21.5 
17.1 
21.6 
17.9 
8.7 
14.0 
10.2 
10.8 
14.9 
8.1 
16.6 
13.8 
18.6 
6.9 
14.5 
19.9 
14.1 
20.4 
14.9 
18.2 
19.4 
19.3 
13.1 
19.5 



ATG20 J023356-503020 

ATG20 J024040-542933 

ATG20 J024313-510517 

ATG20 J024614-495350 
ATG20 J024948-555615 



ATG20 J025329- 

ATG20 J025838- 

ATG20 J030055- 

ATG20 J030328- 
ATG20 J030616- 
ATG20 J031207- 
ATG20 J031425- 



544151 

505204 

510229 

523433 
552808 
554133 
510431 



ATG20 J032327-522630 
ATG20 J032650-533701 
ATG20 J033002-503519 
ATG20 J033114-524148 
ATG20 J033126-525830 



ATG20 J033553-543025 

ATG20 J034349-524116 
ATG20 J034941-540106 
ATG20 J035128-514254 

ATG20 J035700-495547 

ATG20 J040400-552023 

ATG20 J041137-514923 
ATG20 J041247-560035 
ATG20 J041313-533200 

ATG20 J042504-533158 
ATG20 J042842-500534 
ATG20 J042852-543007 
ATG20 J042908-534940 
ATG20 J043221-510925 
ATG20 J043652-521639 

ATG20 J044158-515453 

ATG20 J044748-515100 
ATG20 J044822-504133 
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ACT ID 



RA (J2000) Dec 



S/N 



s a 

(mjy) 



Sdb 

(mjy) 



tint C 

(min) 



AT20G ID 



ACT-S J045103-493632 


01 


51 


03 


-49 


36 


32 


12.1 


58.2 


58.2i« 


ACT-S J045239-530637 


01 


52 


39 


-53 


06 


37 


8.4 


23.3 


22.6tH 


ACT-S J045503-553119 


01 


55 


03 


-55 


31 


19 


9.8 


36.7 


or q + 3.8 


ACT-S J045559-530236 


01 


55 


59 


-53 


02 


36 


6.9 


18.3 


17.4l 2 2 ;? 


ACT-S J050018-532122 


05 


00 


18 


-53 


21 


22 


10.6 


28.2 


27.6±2'7 


ACT-S J050401-502314 


05 


04 


01 


-50 


23 


11 


15.8 


59.6 


59-5±l:l 


ACT-S J051355-505543 


05 


13 


55 


-50 


55 


13 


6.5 


22.0 


20.911;* 


ACT-S J051812-514359 


05 


18 


12 


-51 


43 


59 


8.7 


24.7 


23.9i|8 


ACT-S J052044-550832 


05 


20 


41 


-55 


08 


32 


7.0 


24.1 


22-9±l;| 


ACT-S J052317-530836 


05 


23 


17 


-53 


08 


36 


8.1 


22.9 


22.1+J;? 


ACT-S J052743-542614 


05 


27 


13 


-54 


26 


11 


8.9 


26.0 


25.3i|« 


ACT-S J053117-550431 


05 


31 


17 


-55 


04 


31 


14.1 


47.9 


47 Q+3.3 


ACT-S J053208-531035 


05 


32 


08 


-53 


10 


35 


18.9 


51.2 


5l.2±£;? 


ACT-S J053323-554941 


05 


33 


23 


-55 


49 


11 


16.7 


65.7 


65.71" 


ACT-S J053458-543903 


05 


34 


58 


-54 


39 


03 


10.0 


30.8 


3O.1I31 


ACT-S J053909-551059 


05 


39 


09 


-55 


10 


59 


8.9 


30.6 


29.7l|| 


ACT-S J054025-530354 


05 


40 


25 


-53 


03 


51 


11.5 


31.7 


3i-i±2:i 


ACT-S J054029-535628 


05 


40 


29 


-53 


56 


28 


5.5 


14.8 


13 6+ 2 ' 8 
±o.o_ 2 .8 


ACT-S J054046-541825 


05 


40 


16 


-54 


18 


25 


133.4 


416.2 


416.2 + |-j 


ACT-S J054223-514259 


05 


42 


23 


-51 


42 


59 


27.4 


79.7 


79 7+2-9 
' a, '-2.9 


ACT-S J054357-553206 


05 


43 


57 


-55 


32 


06 


6.2 


22.2 


20.91H 


ACT-S J054830-521836 


05 


48 


30 


-52 


18 


36 


9.3 


25.2 


24.611;? 


ACT-S J054943-524629 


05 


49 


13 


-52 


46 


29 


63.4 


180.3 


I8O.3I28 


ACT-S J055047-530502 


05 


50 


47 


-53 


05 


02 


9.0 


24.5 


23.8l|| 


ACT-S J055152-552642 


05 


51 


52 


-55 


26 


12 


8.6 


31.5 


30.4+3? 


ACT-S J055811-502957 


05 


58 


11 


-50 


29 


57 


16.0 


57.5 


57.511* 


ACT-S J055830-532640 


05 


58 


30 


-53 


26 


10 


6.1 


16.6 


I5.6ll;| 


ACT-S J055946-502656 


05 


59 


16 


-50 


26 


56 


12.5 


45.9 


45.211;? 


ACT-S J060213-542509 


06 


02 


13 


-54 


25 


09 


12.5 


38.3 


37.7i!;J 


ACT-S J060749-525747 


06 


07 


19 


-52 


57 


17 


13.1 


36.0 


OO .0 Q <-r 


ACT-S J060849-545650 


06 


08 


19 


-54 


56 


50 


38.4 


136.9 


136.9±|;t 


ACT-S J061715-530615 


06 


17 


15 


-53 


06 


15 


13.1 


39.6 


39.0+f? 


ACT-S J061846-532948 


06 


18 


16 


-53 


29 
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Flux density as measured directly from the ACT 148 GHz map. 

Deboosted flux densities as described in Section |3. 2. 2l for sources with S m < 50 mjy. For sources measured flux in above 50 mjy, the 
measured flux together with S/N derived errors are reported (See Section l3,2.20 . 



Integration time per square arcminutc 



